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Quantum computation based on quantum cellular automata (QCA) can greatly reduce the control 
and precision necessary for experimental implementations of quantum information processing. A 
QCA system consists of a few species of qubits in which all qubits of a species evolve in parallel. We 
show that, in spite of its inherent constraints, a QCA system can be used to study complex quantum 
dynamics. To this aim, we demonstrate scalable operations on a QCA system that fulfill statistical 
criteria of randomness and explore which criteria of randomness can be fulfilled by operators from 
various QCA architectures. Other means of realizing random operators with only a few independent 
operators are also discussed. 
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The traditional approach to quantum computation has 
been through the circuit model: a series of one and two- 
qubit gates are applied to specified qubits in a specified 
order Such an architecture requires exquisite con- 
trol of each individual qubit and accurate localization 
of the external Hamiltonian. An alternative approach 
to quantum computation utilizes quantum cellular au- 
tomata (QCA). A QCA system consists of just a few 
(typically 1-3) species of qubits such that all qubits of a 
species are addressed equivalently and simultaneously. 

The idea of using a QCA system for quantum compu- 
tation was introduced by Lloyd Q more than a decade 
ago. Lloyd demonstrated the universality of a three- 
species QCA and provided pulse sequences for funda- 
mental gates. Further work on QCA has concentrated 
on proving universality Q, including the universality of 
a two species QCA that is unable to distinguish the left 
neighbor from the right Q . Only recently has there been 
an attempt to exploit the uniqueness of the QCA ar- 
chitecture to enhance quantum information processing 
protocols. Brennen and Williams []| demonstrated the 
production and manipulation of entanglement in a QCA 
architecture. In this work we attempt to utilize the QCA 
architecture in the study of complex quantum dynamics. 

Classical cellular automata (CCA) are systems that 
follow a simple set of local rules applied uniformly to a 
lattice of cells 0. Each cell can have an arbitrary num- 
ber of possible states and, at every time interval, the 
state of each cell is updated based on its current state 
and the state of a given radius of nearest neighbors. For 
example, the evolution of a two-state, radius one, one- 
dimensional CCA updates cell j based on its own state 
and the state of its nearest neighbors. In this case, the 
evolution has eight update directives, one per combina- 
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tion of the three two-state cells, j — 1, j, and j + 1. All 
cells in a CCA evolve in parallel. This is done by copying 
the current CCA state for the cells to refer to when up- 
dating. Though an apparently simple system, CCA rules 
can develop complex behavior and can simulate a wide 
range of phenomena from lattice gases to traffic flow. 

A QCA consists of a lattice of quantum cells, each with 
an arbitrary number of levels. The dynamics of each cell 
can depend on a given radius of nearest neighbors, but 
is restricted by the requirement of unitary dynamics. In 
addition, the no-cloning rule outlaws parallel updating. 
The latter obstacle can be overcome by using at least two 
species of qubits and updating each species separately [j| . 
In this work we assume two-level quantum cells, referred 
to as qubits. In addition, we explore only radius one evo- 
lution, in which each qubit interacts only with its nearest 
neighbors. 

Given that CCA are valuable in the simulation of com- 
plex classical systems, it is natural to ask whether a QCA 
could be used to simulate complex quantum systems. Of 
course, a QCA that is universal can simulate any dynam- 
ics. The question is whether the unique architecture of 
the QCA can provide a less taxing experimental venue or 
provide further insight into complex dynamics. As a first 
step towards answering these questions, we explore the 
ability of a QCA to implement random unitary operators 
efficiently. 

Random matrices were introduced by Wigner to de- 
scribe the energy levels of atomic nuclei 0. Since then, 
random matrices have been used as statistical models for 
a host of complex systems in many areas of physics || . 
From a quantum computation standpoint, some of the 
important systems modeled via random matrices include 
quantum systems whose classical dynamics are chaotic |9|, 
decoherence ^(j, an d quantum computer error models 
Thus, the ability to implement a random unitary op- 
erator allows for the simulation and study of these types 
of systems. 

Beyond simulations, random unitary operators and 
random quantum states, created by applying a random 
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unitary to a computational basis state, play a vital role in 
quantum communication and computation. In quantum 
communication, random states are known to saturate 
the classical communication capacity of a noisy quantum 
channe l 1 121 . In addition, superdense coding of quantum 
states |l3|. a reduction in key length for approximate 
encryption of quantum states, and the construction of 
more efficient data hiding schemes are protocols enabled 
by random unitary operators 0] • Random unitaries can 
also decrease the classical communication cost for remote 
state preparation 0] . Quantum computing protocols fa- 
cilitated by random unitaries include quantum process 
tomography via a fidelity decay experiment using a ran- 
dom operator. Since random unitaries will not commute 
with noise sources effecting the syst em, they can identify 
the type and strength of the noise |l6| . Random quan- 
tum states can be used for unbiased sampling, and the 
amount of multi-partite entanglement in random states 
approaches the maximum at a rate exponential with the 
number of qubits in the system |l7j . 

The appropriate measure against which random uni- 
tary operators and quantum states are defined is the 
Haar measure on the group U(N), where N is the dimen- 
sion of Hilbert space. The random ensemble of unitaries 
drawn from this measure is the circular unitary ensemble 
(CUE) 18]. Unfortunately, an exact parameterization of 
CUE via the Hurwitz decomposition requires expo- 
nential computing resources. Recently, however, pseudo- 
random unitary operators were introduced as efficiently 
implementable substitutes that fulfill various criteria of 
randomness and can be used in the above mentioned pro- 
tocols |l6j. 

The algorithm to produce pseudo-random operators, 
or maps, consists of m iterations of the n qubit gate: 
apply a random SU(2) rotation on each qubit, then evolve 
the system via all nearest neighbor couplings 0. A 
random SU(2) rotation on qubit j during iteration i is 
defined as 
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and the nearest neighbor coupling operator is 

Unnc = e w)Z-°i*°i+\ (2) 

where a 3 z is the z-direction Pauli spin operator. The ran- 
dom rotations are different for each qubit and each iter- 
ation, but the coupling is always 7r/4 to maximize entan- 
glement. After the m iterations, a final set of random 
rotations is applied. 

This paper suggests a modified version of the above 
algorithm that can generate pseudo-random maps appli- 
cable to a QCA. The ability to efficiently generate such 
operators indicates that complex systems can be mod- 
eled and explored on a QCA. Moreover, the number of 



iterations needed to create the pseudo-random operators 
for QCA is comparable the number needed for algorithms 
using a circuit model architecture. 

The modification of the algorithm for application to 
a QCA system is straightforward. For each iteration, 
apply species specific random rotations followed by near- 
est neighbor coupling. Thus, iteration i of a QCA ran- 
dom map consists of applying random rotation Uf on 
all qubits of species A, followed by a different random 
rotation, Uf ' , applied to all qubits of species B, and so 
on for all k species of qubits, followed by U mc , given in 
equation^ In keeping with the original algorithm a final 
rotation of the qubits is always applied. 

This work may be viewed from a different perspective: 
an examination of how difficult (or easy) it is to create 
pseudo-random operators. The algorithm of [lfij requires 
(3mn 4- 1) independent variables for a pseudo-random 
operator. Three independent variables per iteration for 
each random rotation, and one more for the coupling con- 
stants. However, as noted in [lfij . any universal gate set, 
no matter how biased, asymptotically generates the uni- 
form measure of unitary operators. This does not imply 
that a universal gate set will efficiently generate the uni- 
form measure of unitary operators. Nor does it suggest 
that a non-universal gate set cannot display some charac- 
teristics of randomness. Here, we attempt to reduce the 
number of independent variables required for the pseudo- 
random operator algorithm and sec if it is still possible to 
efficiently generate CUE-like statistics. If not, we explore 
whether the generated operators demonstrate any char- 
acteristics of randomness. The modified pseudo-random 
operator algorithm for QCA requires only 3mA; 4- 1 inde- 
pendent variables. Other possible ways of reducing the 
number of independent variables will also be discussed. 

It is important to state that reducing the number of 
independent variable does not necessarily reduce the ex- 
perimental difficulty in realizing the algorithm. Rather, 
experimental realizations would tend not to exhibit the 
symmetries that arise by reducing the number of inde- 
pendent variables as this would require acute precision 
over the internal and external system Hamiltonians. In 
this work, the first attempts at randomness are always 
via operators with maximum symmetries as these are the 
most difficult cases with which to achieve randomness. 

Throughout this paper we avoid specifying the actual 
quantum computing hardware beyond the assertion of a 
a z a J z +1 coupling between nearest neighbor qubits. The 
ai<J z +1 interaction is used in the original pseudo-random 
operator algorithm and is appropriate for certain pro- 
posed realizations of quantum information processing. 
Systems such as quantum dots, however, interact via the 
Heisenberg interaction, s J • s J+1 , which has a total spin 
symmetry. Rotating all the qubits in parallel, as done 
for the k = 1 QCA, commutes with the Heisenberg in- 
teraction and no randomness will be generated. This 
symmetry can be broken and random statistics regained 



3 




FIG. 1: Distribution of the multi-partite entanglement, Q, 
for m = 16 (x), 24, (+), 32 (*), and 40 (o) compared to 
the distribution of Q for CUE random operators (solid line). 
Q is calculated for the wavefunctions produced from evolv- 
ing all (256) possible computational basis states under each 
of 500 eight-qubit maps for each m value. For m = 40 the 
distribution is barely distinguishable from the CUE distribu- 
tion. These distributions are similar to those found in 
despite the further constraints imposed by the QCA architec- 
ture. The inset shows the difference between the average Q 
of states from the QCA operators and operators drawn from 
CUE as a function of m. 



where pj is the reduced density matrix of qubit j. As 
m, the number of iterations, increases, the distribution 
approaches that of CUE operators acting on the same 
states. For m = 40 the distribution of Q is practically 
indistinguishable from that of CUE ^3] • This is the same 
value of m necessary for the circuit model pseudo-random 
operator algorithm to produce a similar distribution. 

Though the maps generated for the QCA, k = 1 archi- 
tecture follow the CUE entanglement distribution, they 
deviate from the expected random statistics for other 
important distributions. Perhaps the most widely used 
statistic for the randomness of operators are the spacings 
between nearest neighbor eigenvalues (or in the case of 
unitary matrices, eigenangles}. For matrices of the CUE 
the expected distribution is 
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where s is the difference between two ordered eigenval- 
ues. The mirror symmetry of the system insures that 
the operator eigenvalues will not follow this distribution. 
Rather, the statistics follow the superposition of two in- 
dependent CUE spectra, as shown in figure @. The total 
distribution for a matrix of two unequal size blocks both 
with CUE distribution, derived from [23J, is 



with a two-species QCA. Moreover, the spin-orbit cou- 
pling present in actual quantum dots, generally consid- 
ered a negative and worthy of cancellation [2(J, introduces 
anisotropy into the system, breaking the symmetry and 
allowing for randomness even with the k — 1 QCA. 

The introduction of different random rotations at each 
iteration is not in accordance with typical CCA evolution, 
which follows the same rule at each time step. Nonethe- 
less, the advantages of the QCA architecture, namely the 
reduced need of external Hamiltonian localization and 
control requirements, are still manifest in this algorithm. 

The first QCA we explore is a single k = 1 species 
chain. Every qubit in the chain rotates in parallel at 
every iteration. This map requires only 3m + 1 inde- 
pendent variables and has an inherent mirror symmetry 
stemming from the equivalence of evolution for all qubits 
barring those at the edges of the chain. Perhaps surpris- 
ing is that, despite the simplicity of the system and the 
inherent mirror symmetry of the operators, the operators 
fulfill several statistical measures of randomness. 

As mentioned above, the production of entanglement is 
one of the motivations for implementing random unitary 
operations. Figure shows the distribution of Q, the 
multi-partite entanglement measure [2lt |22| . for states 
produced from QCA maps operating on computational 
basis states: 
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where gi is the fraction of Hilbert space spanned by block 
% = 1,2, and erf 4 = erf(^tf) and erfc 4 = erfc(^) are 
the error function and complementary error function re- 
spectively. For an 8 qubit operator with mirror symme- 
try, gi = 15/32, (72 = 17/32, and the resulting distri- 
bution is barely distinguishable from the g\ = gi = 1/2 
case, 
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The elements of the eigenvectors of random operators 
also follow ensemble specific distributions. For CUE the 
appropriate distribution, as N — ► 00, is [l8| : 



CUE 
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where y = Nrj, and 77 is the squared modulus of the 
eigenvector element. The randomness of the eigenvector 
elements determines the systems resp onse to perturba- 
tion in the sense of fidelity decay |2J|. The eigenvector 
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FIG. 2: Statistical measures of randomness for 100 one- 
species, m — 40, eight-qubit QCA maps. The left figure shows 
the distribution of s, the nearest neighbor eigenangle spac- 
ings compared to that expected for random unitary operators 
(solid line), regular integrable operators (dotted line), and 
a CUE operator with mirror symmetry (dash-dot line). The 
distribution follows that expected for a random operator with 
mirror symmetry. The right figure shows the distribution of 
the magnitude of the eigenvector elements y compared to the 
distribution expected from CUE. There is a noticeable devia- 
tion from the random distribution. This is seen most clearly 
by the large number of very small terms, y < 10~ 5 , depicted 
by the spike at the left of the figure. Nevertheless, the fidelity 
decay behavior of the maps (not shown), which depends on 
the eigenvector statistics, follows the expected exponential of 
random maps. 



element distribution of k = 1 QCA maps deviate slightly 
from the CUE distribution, as seen in figure @- Nev- 
ertheless, these maps demonstrate the exponential decay 
of fidelity at the rate expected for CUE operators. 

To summarize, while these operators are not random 
with respect to the Haar measure, they appear suffi- 
ciently random for entanglement production. They may 
also be used for protocols relying on randomness of eigen- 
vectors such as random operator quantum process tomog- 
raphy |16| (for noise that does not have the symmetry of 
the maps) and for modeling complex quantum dynamics 
that have inherent symmetries. 

The failure of the above operators to fulfill certain cri- 
teria of randomness can be rectified by breaking the mir- 
ror symmetry of the system. We provide two examples. 
The first is by having a two species QCA chain ABAB . . . 
with an even number of qubits. For this k = 2 map, iter- 
ation i of the circuit consists of a random SU(2) rotation 
Iff on all qubits of species A, followed by a random SU(2) 
rotation Uf applied to all qubits of species B, followed by 
coupling between nearest neighbor qubits, U nnc . Thus, 
the number of independent variables for the generated 
operator is 6m for the two random rotations per itera- 
tion plus one more for the coupling. As shown in figures 
© and au °f the tested criteria of randomness are 



FIG. 3: Distribution of the multi-partite entanglement, Q, 
for two-species QCA maps of m = 16 (x), 24, (+), 32 (*), 
and 40 (o) compared to the distribution for CUE random op- 
erators. As in the case of a one species QCA, the m = 40 
distribution is barely distinguishable from the CUE distribu- 
tion, despite the constraints imposed by the QCA architec- 
ture. The inset shows the difference between the average Q 
for states from the QCA and CUE operators as a function of 



0.8 



0.6 



0.4 



0.2 




0.8 






0.7 






0.6 






§0.5 






§: °' 4 






0.3 






0.2 






0.1 







-4 -2 

log(y) 



FIG. 4: Statistical measures of randomness for 100 two- 
species, m = 40, QCA maps. The left figure shows the distri- 
bution of s the nearest neighbor eigenangle spacings compared 
to that expected for CUE operators (solid line), and regular 
integrable maps (dotted line). The right figure shows the 
distribution of the magnitude of the eigenvector elements y 
compared to that expected for CUE (solid line) . Both criteria 
of randomness are fulfilled by the k = 2 map. 



fulfilled for these maps. 

A second way to break the mirror symmetry of the 
k = 1 QCA is by changing the value of one of the near- 
est neighbor couplings (though not the center coupling). 
In this way the map requires only 3to + 2 independent 
variables while still fulfilling all of the above criteria of 
randomness. An operator generated from a system with 
unequal nearest neighbor couplings is especially relevant 



5 



0.8 



0.6 



0.4 



0.2 




0.8 



0.6 



0.4 



0.2 



-2 

log(y) 




-2 
log(y) 



FIG. 5: Statistical measures of randomness for 100 m = 40, 
one-species QCA rings. In this instance, all nearest neighbor 
couplings are equal to 7r/4 except for two, one of which has a 
coupling of 7r/5 and the other a coupling of tt/4.5. The left fig- 
ure shows the distribution of s the nearest neighbor eigenangle 
spacings compared to that expected for CUE operators (solid 
line) and regular integrable maps (dotted line) . The right fig- 
ure shows the distribution of the magnitude of the eigenvector 
elements y compared to the CUE distribution. QCA architec- 
ture with unequal nearest neighbor couplings (be it a ring or 
chain) are the most likely for experimental implementations. 



for experimental implementations. For many quantum 
computer hardware proposals attaining exact equal cou- 
plings between qubits is nearly impossible. An actual 
chain of qubits, such as in nuclear magnetic resonance 
or quantum dots, could not be expected to have the mir- 
ror symmetry assumed above. As we have demonstrated, 
however, this allows such systems to more easily generate 
pseudo-random states and operators. 

The discussion to this point has centered around QCA 
chains. If the QCA were structured as a ring (which, ex- 
perimentally, may be more difficult) such that all qubits 
had two nearest neighbors, we have checked that the sym- 
metries inherent in the system do not allow for enough 
entanglement production following the random distribu- 
tion of Q, regardless of the number of iterations. How- 
ever, if one of the couplings is (even slightly) different 
than the others the random distribution of Q is recov- 
ered as the system is then equivalent to the QCA chain. 
If two couplings are different (from the others and each 
other) , all symmetries have been broken and full pseudo- 
randomness is recovered as shown in figure JSJ. 

There are other ways to reduce the number of indepen- 
dent variables in the pseudo-random operator algorithm 
besides a QCA. We explore some of these possibilities in 
attempt to achieve aspects of randomness with as few in- 
dependent variables as possible. First, we explore what 
is, in some sense, the opposite of the QCA discussed 
above. For the QCA operators the same rotation was 
applied to every qubit but the rotation was different for 



FIG. 6: Statistical measures of randomness for 100 m = 
40, repeat maps. The left figure shows the distribution of 
s the nearest neighbor eigenangle spacings compared to the 
expected distributions for CUE operators (solid line), regular 
integrable systems ( dotted line) and COE operators (dashed 
line) . The right figure shows the distribution of the magnitude 
of the eigenvector elements y compared to the CUE (solid line) 
and COE (dashed line) distributions. Both distributions are 
similar to those of the COE universal symmetry class. 



each iteration. We now study operators in which a dif- 
ferent random rotation is applied to each qubit but the 
rotation is the same at every iteration. As before there 
is a 7r/4 nearest neighbor coupling evolution between the 
rotations. These operators, to which we shall refer to as 
repeat maps, require only 3n + 1 independent variables. 

Given the structure of repeat maps, the repetition of 
the same operation over and over, it is clear that they 
have time-reversal invariance, applying the map back- 
wards produces the same results as applying it forward. 
Hence, the eigenvalue and eigenvector statistics deviate 
only slightly from the distribution appropriate for ran- 
dom orthogonal matrices. The circular orthogonal en- 
semble (COE) includes unitary operators that have anti- 
unitary symmetry (time reversal invariance) and is a sub- 
set of t he g eneral CUE. The level spacing for the COE 
class is iTSl: 
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and the distribution of the elements of the eigenvectors 
of COE matrices as N — ► oo are (lsf 
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The operators of many quantum analogs of classically 
chaotic systems belong to the COE class and, therefore, 
repeat maps may form appropriate models for these sys- 
tems. However, repeat maps do not produce the entan- 
glement distribution expected for random operators. 

Finally, we turn to versions of the pseudo-random op- 
erator algorithm which faithfully follow the evolution of 
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CCA. As explained, there are two major characteristics 
of CCA evolution: homogeneity of evolution for each cell, 
and homogeneity of evolution at each time step. In the 
first part of this work, we examined evolution in accor- 
dance with only the first of these characteristics. Repeat 
maps describe evolution with only the second character- 
istic. Currently, we explore the dynamics of maps that 
evolve as CCA in both respects: each qubit rotates via 
the same random SU(2) rotation, U, and that rotation 
is the same for all iterations. This gives a total of only 
3+1 independent variable for the entire operator. 

As these QCA maps are even more limited than the 
repeat maps discussed above, it is no surprise that the 
entanglement produced by these maps do not follow the 
distribution of random maps. However, these maps de- 
viate only slightly in the other criteria of randomness. 
As shown in figure (JJJ the eigenvector element distribu- 
tion deviates somewhat from the COE distribution while 
the nearest-neighbor spacing distribution deviates only 
slightly from the global statistics expected from a map 
with two differently sized COE blocks due to mirror sym- 
metry [3 
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with <?i and gi defined as above. As with the CUE mirror 
symmetry operator, the resulting distribution is barely 
distinguishable from the g\ = g^ = 1/2 case 
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As with the QCA maps discussed at the beginning of 
this work, the mirror symmetry of the system can be 
broken by changing one of the couplings. In this way the 
eigenvalue and eigenvector statistics revert back to the 
COE distributions seen in the repeat map and shown in 
figure © • This operator requires only 3 + 2 independent 
variables. 

Perhaps it should come as no surprise that operators 
with so few independent variables can still fulfill crite- 
ria of randomness. Quantum chaos models, such as the 
quantum sawtooth and quantum Harper's map [2H Eif , 
have only one or two free parameters, yet, fulfill criteria 
of randomness, and can be efficiently implemented on a 
quantum computer. What we have shown here, however, 
is that quantum chaos models are not exceptional cases, 
showing randomness due to their connection with a clas- 
sically chaotic analog. Rather, most operators with few 
random variables will still show many characteristics of 
randomness. Regularity is the exception, randomness is 
the general rule. 
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FIG. 7: Statistical measures of randomness for 100 m = 40, 
k — 1, QCA maps in which there is only one random rotation. 
This one rotation is applied to each qubit at each iteration. 
The left figure shows the distribution of s the nearest neighbor 
eigenangle spacings compared to the global spectral statistics 
of a COE operator with mirror symmetry (dash-dot line), 
and the right figure shows the distribution of the magnitude 
of the eigenvector elements y. Both distributions are similar 
to those of the COE universal symmetry class. 
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FIG. 8: Statistical measures of randomness for 100 m = 40, 
k — 1, QCA maps in which the same random rotation is 
applied to each qubit at each iteration. In these maps one of 
the nearest neighbor couplings in the chain is tt/5 while all 
the rest are 7r/4. The left figure shows the distribution of s 
the nearest neighbor eigenangle spacings compared to COE, 
and the right figure shows the distribution of the magnitude 
of the eigenvector elements y compared to that expected to 
COE. 



In conclusion, we have demonstrated that general as- 
pects of complex quantum dynamics can be studied on 
a QCA architecture. The evidence behind this supposi- 
tion is the ability to generate pseudo-random operators 
that fulfill general criteria of randomness. These oper- 
ators are appropriate substitutes for random operators 
which are often used as models of complex quantum sys- 
tems and play an important role in other quantum infor- 
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mation processing protocols. Throughout this work we 
have attempted to introduce operators generated by as 
few independent variables as possible. These operators 
are tested to determine whether they meet the various 
criteria of randomness. Operators which do not fulfill all 
of the criteria may nevertheless prove useful in certain 
simulations of quantum systems and other quantum com- 
putational protocols. The minimum number of indepen- 
dent variables is reached via simulation of an algorithm 
which mimics classical cellular automata evolution. Yet, 
even these operators demonstrate many aspects of ran- 
domness. This implies that even with few independent 
variables most operators will tend towards randomness. 
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